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Abstract 

The extraction of the CKM angle a from the asymmetry in B° — > 7r + 7r _ vs -B — > 
7r + 7r~ suffers from a currently unknown penguin contribution. Experimentally, one 
can determine the magnitude and phase of the CP asymmetry from a time-dependent 
analysis of tagged events, and the average rate for B° and B° decays to it + -k~ from 
untagged events. These measurements, together with the magnitudes and relative 
phase of the tree and penguin diagrams, can in principle completely determine a, free 
of discrete ambiguities. We perform an error analysis on a given assumptions on the 
values and uncertainties of both the measurements and theoretical inputs. 
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1 Introduction 




The unitarity of the CKM matrix corresponds to six independent conditions between 
the CKM matrix elements 0. Geometrically, 
these conditions can be visualized as triangles 
formed by the appropriate CKM matrix ele- 
ment combinations. The condition on the d 
and b rows, J2 q = u ,c,t VqdVqb = 0; results in the 
CKM triangle shown in Fig. 1. If CP viola- 
tion occurs via the CKM matrix, this triangle 
will have non-zero area. The three angles in 
this triangle, known as a, (3 and 7, can be ex- 
tracted via CP asymmetries: the angle (3 from 
the decay B° — > J/ipK s ^, the angle a from 
the decay B° — > 7r + 7r~[|3| and the angle 7 from 
the decays B s -> D 8 K$^ and B ± -> D Q K±^. 

The amplitude for B — > nn is dominated by the tree process b — > uW, W —>■ du, which 
has a weak phase 7. If this were the only contribution, the CP asymmetry in B° — > 7r + 7r~ 
would cleanly measure sin(2/3 + 27) = — sin(2a). However, there may be significant contri- 
butions from the penguin process b — > dg. This process has a weak phase —(3 and therefore 
results in a distortion of the CP asymmetry. Any interpretation of the CP asymmetry in 
B — > 7r + 7r - must therefore account for the possibility of a penguin contribution. 

Gronau and London || have shown that the penguin contributions can be isolated by 
applying an isospin analysis to the decays B° — > 7r + 7r~, B + — > 7r + 7r° and B° — > 7r°7r°. 
Aleksan, Gaidot, and VaisseurQ have estimated that this analysis typically results in a 60% 
increase in the uncertainty on sin(2a) relative to the ideal case where only the tree diagram 
need be considered. While feasible for an experiment at an e + e~ collider, this analysis is of 
no use to an experiment at a hadron collider, given that it is very unlikely that the mode 
B° 



v cd v; b 

Figure 1: The d,b CKM triangle 



7T 7r° will ever be reconstructed in such an environment. 



Silva and Wolfenstein have shown 0] that a can be determined from the CP asymmetry 
in B° — > 7r + 7r~, the relative rates for B° — ► 7r + 7r~ and B° — > K + n~, and assuming SU(3) 
symmetry and factorization. Some complications are the possibility of final state phase 
shifts, and electroweak penguins that invalidate the SU(3) correspondence ||. 

We present herein an analysis of the expected uncertainty on a, and the number of 
discrete solutions, given a measurement of the time-dependent asymmetry between B° — > 
7r + 7r~ and B° — > 7r + 7r~, a measurement of the average branching ratio for B° and B° decays 
to 7r + 7T - , and constraints on the magnitudes and relative phase of the tree and penguin 
diagrams. 
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2 CP violation in B° tt+tt" 

The mathematical expression of the CP asymmetry in the decay B° — > 7r + 7r~ can be found 
in numerous places in the literature. Here we follow the exposition of Gronau||. 

In general, the two physical states Bl and Bh are given in terms of the strong eigen- 
states B° and B° via 

\B L ) = p\B°} + q\B } (2.1) 
\B H )=p\B°)-q\B°) (2.2) 

If two amplitudes (e.g. tree-level and penguin) contribute to the decay B° — > f, then the 
decay amplitudes of B° and B° to a CP eigenstate / are given by 

A f = A{B° -> /) = A T e i<t>T + A P e^ p (2.3) 
A/ = A(B° -»• /) = A T e"^ T + A P e"^ p (2.4) 

where each term in the above expression corresponds to a process. The amplitudes Ai are 
complex and contain hadronic final-state-interaction phases. The time-evolution of states 
initially pure in B° and B° are then given by 

r(£°->/) = IA/IV* [|A| 2 sin 2 (xt/2) + cos 2 (xt/2) -7mA sin art] (2.5) 
T(B° -> /) = I^IV* [sin 2 (a;t/2) + |A| 2 cos 2 (xt/2) + /mAsinxtj (2.6) 



where 



A = ^ (2.7) 



The time-dependent asymmetry, a(t), is thus 

_ r^ -> /) -r(B° -> /) _ (1 - |A| 2 )cosxt-2JmAsinxt 

a(t) ~ r(Bo^/) + r(B°^/) ~ TTpf ( } 

For the case / = , q/p = e~ 2tf3 and 0r = 7, arid in the approximation of neglecting the 

penguin contribution, i.e. Ap = 0, A is a pure phase which results in Im\ = — sin(2/5 + 27) = 
sin(2a), assuming the unitarity of the CKM matrix. In this case, the amplitude of the 
asymmetry directly yields a clean extraction of the angle 2a - but with a discrete ambiguity. 

In this same decay mode, however, the penguin, assumed to be dominated by the top 
quark loop, has a CKM phase given by <pp = —f3 and therefore the extraction of 2(j3 + 7) 
is not clean. Inspection of equation [2.8| shows that the effect of the penguin contribution is 
the addition of an additional sinusoidal modulation in the time-dependent asymmetry, the 
additional factor (1 — |A| 2 ) cosxt. The overall asymmetry can then be written as 

a(t) = Asin(xt + <p) (2.9) 
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where 



1 - |A| 2 ) 2 + 4/m 2 A 
A = 1+ x sign(-imA) (2.10) 



1-|A| 2 7T 



7T 



tm * = ^r--2 < * < 2 (2 - n) 

This convention reduces smoothly to the standard expression for the no-penguin case. We 
also exploit the dependence of the average branching ratio for B — > tt + tt~ , B avg on the angle 
a: 

B avg OC T(B° -> TT+tT) + T(B° -> TT+TT) = (A/1 2 [l + |A| 2 ] (2.12) 

In what follows we will refer to the strength of the penguin contribution, Ap, relative 
to the tree-level, A T . We thus introduce the ratio of the amplitudes / = A P /A T = \f\e l5 , 
where 5 is the strong phase difference between the amplitudes, and obtain for A: 

1 — fp~ ia 

A = e 2ia (2.13) 

1 - fe ia v ; 

Experimentally, we have three observables: the magnitude of the asymmetry, A, the 
phase of the asymmetry at t — 0, (ft, and the average branching ratio, B avg . All of these 
are affected by both the tree-level and penguin diagrams. The dependence of the three 
observables, A, <fi, B avg , on the angle a is shown in figure [|, for various values of 5. We have 
taken |/| = 0.2 for these plots. We see that the asymmetry in the presence of a penguin 
contribution is no longer symmetric around a = 90°. Also, at a = 90°, the average branching 
ratio is no longer equal to the tree-level one, but it is increased by a factor 1 + \ f\ 2 , i.e. 1.04 
in this example. Note that the phase, 0, vanishes for 5 = 0, 180°. 

In the absence of penguins, there is a four-fold ambiguity in the determination of a from 
the CP asymmetry: The asymmetry is identical for the case a — > n/2 — a, and for a — > a — i\. 
A most interesting feature of the plots in figure || is the behavior of the branching ratio 
between and 7r/2: the curves change monotonically and thus lift the ambiguity between a 
and 7r/2 — a. Also, the curve for A is antisymmetric around a = 0, while the curves for and 
B avg are symmetric around a = 0. Thus, the ambiguity between a and a — 7r is also lifted: 
there is only a single discrete case, A = and a —>■ —a, where two values of a are a solution 
given A, <fr, and B avg . In summary, in the presence of penguins, the four-fold ambiguity is 
in principle completely lifted except for the discrete case A = where it becomes a two-fold 
ambiguity, a — > — a. 

In the next section we estimate the expected error on from fitting the above asym- 
metry as a function of the statistics. Most proposals for experiments at hadron colliders 
involve a 7r + 7r~ trigger that imposes (usually indirectly, via impact parameter requirements) 
an effective cut t > T on the lifetime of the B° decays reconstructed. We thus compute the 
error on the observables as a function of the effective cut value, T. 
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Figure 2: The three experimental quantities, the asymmetry, A, the phase, 0, and the average 
branching ratio, B avg (normalized to the "no penguin" case), as function of the angle a. The 
relative size of the penguin contribution is 20%. The various curves correspond to different 
values of the strong phase difference 5. 
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3 Fitting the data for A and <j) 



The numbers of B° (N+(t)) and B° (iV_(t)) at time t can be written as 

N±(t) = — e~* [1 ± A sin(xt + 0)] (3.14) 



where we used equation |2l| and integrated over all time to express \A\ 2 in terms of the total 
number of B° and B° mesons, N. We are interested in estimating the error on the quantities 
A and 4> resulting from a fit to the data by the above form. The probability that a set of B° 
and B° mesons (initially pure) will decay at times ti and tj respectively is given by 

C = IT e ~ U l 1 + A Sin ( xt i + 0)] IT e_<J t 1 - A Sin ( xt 3 + 0)] ( 3 - 15 ) 
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Figure 3: Contours of equal probability in the^A — 5(j) plane 
Equation |3.15| is then the expression for the likelihood for this set of events. The inverse 



of the covariance matrix for the variables A and 0, G = V l , is given by 



G 



AA 



d 2 ln£ 

dA 2 



a 2 ln£ 



G 



Ad> 



d 2 ln£ 



(3.16) 
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The covariance matrix is computed in appendix A. 

The one sigma contour in the A — plane is given by the ellipse equation: 



1 



P 1 



5A 2 2p5A5<p 



+ 



(75 



(3.17) 



where 5 A = A — A, 5<p = <f) — <f) and A, <j> are the "true" values, and p is the correlation 
coefficient. This error ellipse is shown in figure |3] for four different values of the asymmetry 
A, and = 0.1. It can be seen that the error on <p decreases as the asymmetry gets larger. 



4 Extraction of a given A, 0, and B avg 

As discussed in Section |2|, there are three observables related to the CP asymmetry in B — ► 
7T7T. The uncertainties on the measurements of A and were discussed in Section [| The third 
observable, B avg , can be measured from untagged events, and will therefore have a very low 
statistical error. The extent to which systematic uncertainties can be controlled will therefore 
be a crucial consideration. While absolute branching ratios are very difficult to obtain, it 
will be sufficient to measure the branching ratio relative to the process B° — > £ + tt~v. 
There are five unknowns that influence the values of the 3 observables: 

• a, the weak phase we are trying to extract from these measurements. 

• At, the amplitude of the tree diagram. Assuming factorization, the amplitudes for 
B° — > l + n~v and B° — > 7r~7r + are proportional to a common form-factor, spanning a 
range of q 2 for the first case, and evaluated at q 2 = m 2 ^ for the latter casefllQfl . Ref. |~l| 
points out that color-allowed B decays are well described by factorized amplitudes, but 
find that it is necessary to add non-factorized amplitudes to describe color-suppressed 
B decays. Since the decay B° — > 7r~7r + is color allowed, we assume that the decay 
B° — > £ + 7T~b> will be observed in conjunction with B° — > 7r~7r + and used to predict 
A T . 

• Ap, the amplitude of the penguin diagram. This amplitude can be estimated from 
measurements of the decay B° — > K~tt + , applying SU(3) corrections, and scaling by 
|Kd/Vt s |@. Some complications are that this decay may in turn have a contribution 
from tree diagrams, and furthermore, there may be electroweak contributions that 
invalidate the SU(3) correspondence |J. The decay mode B s — > <pp° may possibly be 
used to check our understanding of these effects |]12|- We assume that the relative size 
of the signals for B° — > K~tt + and B° — > 7r~7r + will determine Ap, although with less 
precision than for At- 

• The weak phase of the penguin amplitude. The top quark dominates in the loop, 
therefore this phase will be —f3 to a good approximation! 13| . As shown in equation [2.13| , 
in this case we are not sensitive to the value of (3. 
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• The strong phase difference, 5, between the tree and penguin diagrams. There is a 
perturbative phase difference of order 10° ||, and there are also nonperturbative effects 
from hadronization that are expected to be small but are incalculable. The phase 
shift between the 1=0 and 1=2 amplitudes can be obtained from a measurement of the 
branching ratios of B — > 7r + 7r~, B — ► 7r + 7r°, and B — > This check can be done at 

a symmetric e + e~ collider, and does not require flavor tagging or time-ordering. The 
question then becomes: Is there a difference in the hadronization for a penguin diagram 
and 1=0 tree diagram? The extent to which these phase shifts can be constrained helps 
constrain our solution for a. 



Given an assumption for the central values of these parameters, we can calculate the 
values of A, (ft, and B avg . Given an assumption for the effective number of tagged events, 
Ntaa fll5 1 , we can calculate the error matrix for A and <fi. We will make assumptions on the 
uncertainties on At, Ap, 5, and B avg , parametrized as Gaussians with widths a (At), a(Ap), 
cr(5), and a(B avg ). With these assumptions, we can form a x 2 ■ The minimization program 
MINUIT is used to minimize this x 2 ■, return the input value of a, and estimate the expected 
uncertainty on a measurement of a. 

Unless specfied otherwise, we use the following as default values of the parameters: 



• N tag = 100, starting at cr = 1.6 lifetimes (T^j. 

• A T = 1.0 

• A P = 0.2 

• 5 = 0.0 

In Fig. ^ we show the expected la errors as a function of a for various conditions. In 
Fig. ||a, we show the errors for the case where Ap = 0.0, and has been constrained to zero in 
the fit. We see that the errors are largest for a near 45° and 135°, where the dependence of 
sin(2a) on a is lowest. In Fig. ||b, c, and d, we show the errors for the case where Ap = 0.2 
and we are able to put the specified constraints on the amplitudes. We see that in the case 
where there is a penguin amplitude, and it is well understood, in general, the errors are 
smaller than in the case of no penguin amplitude. 

The plots in Fig. [| do not convey all the relevant information on the constraints on a. 
The errors are highly non-Gaussian, and there are multiple minima. To gain further insight, 
we choose two particular input values for a, 47° and 67°. We then function of the 

assumed value of a in the fit. For each point in the scan, we hold a fixed, and mimimize 



the x 2 with respect to all the other parameters. We then plot J x 2 — Xmin as a function of 



a, interpreting \jx 2 ~ Xmin as t ne number of standard deviations on a. 

We show the results in Fig. || for the case Ap/At = 0.05. In Fig. |]a, we assume 
very little knowledge of the tree and penguin parameters. We see that when the penguin 
contribution is small, loose constraints are sufficient for the determination of a. However, 
even with the tight constraints of Fig. He we are unable to lift the discrete ambiguities. 
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Figure 4: One a errors on a as a function of the input value of a. The solid line shows the 
positive errors, and the dashed line the negative errors: 

(a) Ap = 0, and fixed to in the fits. 
The following assume Ap = 0.2: 

(b) a{Ap)/A P = 0.1, a(A T )/A T = 0.03, a(S) = 20°, and a(B avg )/B avg = 0.06. 

(c) A P , A T , and 5 held fixed in the fits, and a(B avg ) j B avg = 0.06. 

(d) a{Ap)/A P = 0.1, a(A T )/A T = 0.015, a (5) = 20°, and a{B avg )/B avg = 0.03. 
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Figure 5: Na as a function of a, for an input value of a = 47° on the left and a = 67° on 
the right. The input value for Ap/A? is 0.05. 
For the following, we assume a(B avg ) / B avg = 0.06: 

(a) a(Ap)/A P = 1.0, a(A T )/A T = 1.0, no constraint on 5. 

(b) a(Ap)/A P = 1.0, a(A T )/A T = 0.03, no constraint on 5. 

(c) a(Ap)/Ap = 1.0, a{A T )/A T = 0.03, a(S) = 20°. 

(d) a(Ap)/A P = 0.1, a(A T )/A T = 0.03, no constraint on 5. 

(e) a(Ap)/A P = 0.1, a(A T )/A T = 0.03, <t(5) = 20°. 
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Figure 6: Na as a function of a, for an input value of a = 47° on the left and a = 67° on 

the right. The input value for Ap/A? is 0.2. 

For the following, we assume a(B avg ) / B avg = 0.06: 

(a) a(Ap)/A P = 1.0, cr(A T )/A T = 1.0, no constraint on S. 

(b) a(Ap)/A P = 1.0, a(A T )/A T = 0.03, no constraint on 5. 

(c) a(Ap)/A P = 1.0, a(A T )/A T = 0.03, a(S) = 20°. 

(d) a(Ap)/A P = 0.1, a(A T )/A T = 0.03, no constraint on 5. 

(e) cr(Ap)/Ap = 0.1, a{A T )/A T = 0.03, <r((J) = 20°. 
For the following, we assume cr(B avg )/B avg = 0.03: 

(f) a(Ap)/A P = 0.1, a(A T )/A T = 0.015, a(5) = 20° 
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We show the results in Fig. |5| for the case Ap/A? = 0.2. In Fig. |G]a, we assume very 



little knowledge of the tree and penguin parameters. As qualitatively pointed out in Ref [14 



in this case, we can rule out only a small fraction of the available parameter space. As we 
add constraints in Fig. |6|b, c, and d, there are fewer minima and more of the parameter space 
can be ruled out. As shown in Fig. ||e, it is not until we constrain all the parameters that we 
are left with a single minimum. Tightening the constraints in Fig. ^f, we more convincingly 
rule out alternative minima and improve the precision of the measurement of a. 

In Fig. [7], we show similar plots as for Fig. |6|, using the same input values for a but 
negative assumed values for a. As discussed earlier, there is another solution only for the 
discrete case A = 0. However, even for other values of A, two solutions may be allowed within 
the uncertainties of the measurement. Fig. |7] illustrates how well we can choose between the 
solutions. The separation is quite convincing for a = 47°, i.e. for large values of A, and 
becomes more difficult for a = 67°, i.e. as A gets smaller. 



Input a = 47° Input a = 67° 
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Figure 7: Na as a function of —a, for an input value of a = 47° on the left and a = 67° on 
the right. For the following, we assume a(B avg ) /B avg = 0.06: 
(a) a(Ap)/A P = 0.1, a(A T )/A T = 0.03, a{5) = 20°. 

For the following, we assume a(B avg )/B avg = 0.03: (b)a(A P )/A P = 0.1, a(A T )/A T = 0.015, 
a (5) = 20°. 

It is clear from the above plots that knowledge of the strong phase difference 5 decreases 
the error on a and in addition can help distinguish between the discrete solutions on a. 
Figures |] and | are similar to Figure ||, except we scan as a function of assumed values for 

both a and 5, and plot constant contours in \Jx 2 ~ Xmin- m ^ ne case of a = 47° there are 
only two minima, and one needs to know 5 to better than 60° in order to distinguish between 
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Figure 8: One, two, and three a contours as a function of a and 5, for input values of 
a = 47°, S = 0°, and A P /A T = 0.2. 
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Figure 9: One, two, and three a contours as a function of a and 5, for input values of 
a = 67°, 5 = 0°, and A P /A T = 0.2. 
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the two discrete solutions by 2a. In the more complicated case of a = 67°, there are multiple 
solutions even within the first two quadrants for a. Here a more precise knowledge of 5 is 
required to unambiguously extract a. 



a((5) = GO a(<5) = 40° a(<5) = 20 




100 100 100 



Input a, degrees 



Figure 10: One, two, and three a contours as a function of input value of a and assumed 
value of a, for various values of ag. Note that there always remains a discrete ambiguity 
between a and —a where A = 0. 

The above two values of a of 47° and 67° lead to quite different conclusions on the 
knowledge of 5 required for a unique determination of a. We have also investigated the effect 



of the uncertainty on S on the extraction of a, for all values of a. As an example, Figure 10 



shows constant contours of \jx 2 ~ Xmin i* 1 the plane of input value of a and assumed value 
of a, for three different constraints on 5. We see that the discrimination between discrete 
ambiguities improves with the precision on 5 until the uncertainty on 5 reaches a value of 
20°. We find that this result holds for all values of a. 
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5 Conclusions 



We have presented the results of an error analysis on a, given a measurement of the time- 
dependent asymmetry between B° — > tt + 7T~ and B° — > 7r + 7r~, a measurement of the average 
branching ratio for B and B decays to 7r + 7r~, and constraints on the magnitudes At and Ap 
and relative phase S of the tree and penguin diagrams. While there are an infinite number of 
possible scenarios, our results set the scale for how well these parameters need to be known. 
We have considered scenarios in which we have an effective number of 100 perfectly tagged 



signal events) 15 



In the case where the penguin contribution is small (Ap/Ap ~ 0.05), only crude in- 
formation on and At is needed for the extraction of a. However, we are left with four 
discrete ambiguities, as in the case where there is no penguin contribution. 

In the case where the penguin contribution is larger (Ap/Ap ~ 0.2), more precise 
information on and At is needed. If this precision can be achieved, the uncertainty 
on a is in many cases smaller than for the case of no penguin amplitude. Furthermore, 
if it is possible to place constraints on 5, some or all of the discrete ambiguities may be 
lifted. A large penguin amplitude therefore presents an opportunity for a much improved 
determination of a. 

In summary, we find that if the penguin amplitudes are either small or well understood 
then it is possible to determine a from the CP asymmetry in B° — > 7r + 7r~ without resorting 
to the observation of final states with neutral particles. Thus, measurements would be 
feasible at hadron colliders as well as e + e~ colliders. Furthermore, a large penguin amplitude 
presents an opportunity for improved precision on a while lifting some or all of the discrete 
ambiguities. 
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A Likelihood formalism for the extraction of fitting er- 
rors 

The analysis follows closely that in reference [I7|. For brevity, we define the functions n±(t): 

n±(A,(j),t) = 1 ± Asin(xt + 0) (A.l) 

With this notation, and ignoring an irrelevant constant term, the log-likelihood (see equa- 
tion |3.15|) is 



In £(A, 0) = ^ ln n + (A <t>, U) + ln n - ( A > & tj) 



(A.2) 



The second derivatives of the above function are the elements of the inverse of the correlation 
matrix, G = V" 1 . For example, 



d 2 ln C ^ sin 2 (xtj + 4>) ^ sin 2 (xt,- + (f>) 
Gaa = tttt- = 2^ u 7Tw + <L 



These sums can be approximated by integrals: 

E/±(**)» / U{t)N±{t)dt 



(A.3) 



(A.4) 



where the limits of integration assume that a lifetime cut is imposed on the reconstructed 
mesons, and N±{t) are defined in equation |3.14 Since in general f±{t) = g±(t)/[n±(ti)} 2 , 
equation [A.3| (for the general error on variables p and q) thus becomes 



G 



pq 



9+V) + 9-(t) 



n + (t) Ti-(t) 



dt 



2 Jt 

N f°° 

-J t e- t [^ + (t)n_(t)+^(t)n + (*)]dt 



where we have ignored terms which of order A 2 : 

n + (t)n_(t) = 1- A 2 sin 2 (xt + 



(A.5) 
(A.6) 

(A.7) 



We expect this to be a reasonable approximation even for extreme values of sin(2a), since 
in practice the observed asymmetry will be reduced by a dilution factor D ~ 0.5 resulting 
from imperfect flavor tagging. With some algebra, we obtain 



POO 

G aa = N e~ l sin 2 (xt + <j>)dt 

JT 

roc 

/ e~* cos 2 (xt + 4>)dt 

JT 

A e~ t sin2(xt + <f>)dt 

JT 



G^ 
Ga4> 



NA 2 

N 
~2~ 



(A.8) 
(A.9) 
(A.10) 
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We note that, as expected, these equations are invariant with respect to the following trans- 
formation in the presence of dilutions: 



A — > DA (A.ll) 
N 

N - W2 (A.12) 
For brevity, we introduce two new functions, P c (x,T) and P s (x,T) given by 

P c (x,T) = — l — [cosxT - xsinxT] (A. 13) 

1 + x 2 

P a (x,T)= 1 2 [x cos xT + sin xT] (A. 14) 

J- ~~ \~ Ob 



and the end result is 



Gaa = ye" T [1 - P c (2:r,t )] (A. 15) 

G H = A 2 ^-e- T [l + P c (2x,t )} (A.16) 
G A4 > = je- T AP s (2x,t ) (A. 17) 



where t Q = T + 
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